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ABSTRACT 


The  ability  to  measure  and  predict  atmospheric 
turbulence  affecting  laser  beam  propagation  is  a  major 
concern  when  considering  military  applications.  Such  a 
method  using  a  telescope,  high  resolution  television  camera, 
video  tape  recorder,  digital  storage  oscilloscope,  and 
calculator  system  has  been  devised,  tested  and  utilized.  A 
laser  beam  signal  is  recorded  on  video  tape  for  further 
processing.  This  signal  is  displayed,  stored  and  digitized 
using  a  Tektronix  4S3  digital  storage  oscilloscope.  The 
digitized  signal  is  sent  to  a  Hewlett-Packard  9825  computing 
system  for  Fourier  transform  analysis  and  determination  of 
the  refractive  index  structure  constant,  C n2 •  Several 
trials  were  conducted  using  He-He  and  Ga-As  lasers.  The 
results  demonstrated  good  correlation  with  theoretical 


predictions  as  well  as  previously  analyzed  data 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  increasing  use  of  lasers  and  laser  technology  for 
military  applications  has  brought  about  a  need  for  analysis 
of  the  laser  beam  in  its  environment,  the  turbulent 
atmosphere.  A  pro jar t  at  Naval  Postgraduate  School  dealing 
with  this  subject  is  continuing  and  is  the  main  topic  of 
this  thesis. 

The  patterns  produced  by  lasers  on  targets  have  inherent 
problems  that  includa  broadening,  beam  wander  and  intensity 
fluctuations  brought  about  by  turbulence  in  the  atmosphsr®. 
These  effects  of  atmospheric  turbulence  on  laser  propagation 
have  been  well  deterained  [Ref.  1].  In  terms  of  the  Fried 
model,  [Ref.  2],  Cn2,  the  refractive  index  structure 
constant,  has  been  determined  to  adequately  express  the 
Modulation  Transfer  Function  (MTF)  ,  or  the  Mutual  Coherence 
Function  (MCF)  for  the  atmosphere. 

A  system  has  been  developed  that  provides  a  measurement 
of  c  n  2  for  atmospheric  turbulence  along  the  optical  path 
through  which  a  laser  is  propagated.  The  system  employs  a 


vidicon  and  telescops  as  the  detector  and  a  distant  laser  as 
the  source.  This  system  duplicates  the  slit  scancinq  system 
presently  in  use  at  Naval  Postgraduate  School  [Hef.  2]. 
Tests  using  the  vidicon  equipment  have  been  previously 
completed  with  values  for  Cn2  on  the  order  of  10"15  m-2/3 
being  obtained  [Ref.  3].  Measurements  have  been  made  using 
the  same  experimental  set-up,  except  that  the  Tektronix  468 
Digital  Storage  Oscilloscope  are  used  for  digitization 
instead  of  the  Quantex  DS-30  Digital  Video  Processor.  These 
measurements  have  demonstrated  values  of  3n2  of  comparable 
ac curacy. 

B.  PROBLEM 

By  using  the  digital  storage  capability  of  the  Tektronix 
468,  the  previously  recorded  Laser  signal  is  used  as  an 
input  and  is  evaluated  by  modifying  the  program  developed  by 
Crager  [Ref.  4].  \  brief  overview  of  the  procedure  is 

described  below.  Detailed  explantions  of  the  experiment  and 
computer  program  are  contained  in  Chapters  III  and  IV, 
respectively. 

The  approach  taken  is  basically  the  same  as  that  of 
Crager,  but  because  of  the  unavailability  of  a  disk  RDM  it 
is  necessary  to  store  both  raw  and  processed  data  cn 


magnetic  taps.  me  ait  a 


necessary  by  the  H?  93  25, 


_oaaea  and  recorder  as 


The  basic  assumption  permitting 


analysis  is  that  the  horizontal  TV  scan  line  through  the 
laser  spot  is  considered  to  accurately  mirror  a  point  spread 
function  of  the  image.  Pinal  analysis  using  tha  468 
oscilloscope  has  shown  this  to  be  a  valid  hypothesis.  The 
digitized  data  from  the  468  oscilloscope  agree  with  the 
previously  measured  lata  which  use  many  pixels  of  television 
data  digitized  by  the  DS-30  to  express  the  point  spread 
function. 

The  sequence  of  analysis  is  that  a  video  tape  recording 
is  made  of  the  TV  image  of  the  propagated  laser  beam.  The 
output  of  tha  video  recorder  is  sent  to  the  Tektronix  468, 
where  the  derived  TV  scan  line  is  digitized,  averaged,  and 
stored.  The  HP  9825  records  the  ligitized  data  and  produces 
a  line  spread  function  (LSF)  ,  by  integrating  the  point 
spread  function. 


y  (max) 


L3  F  (x)  =  |  PSF  (r)  iy 


(1.1) 


=  V<x2+7a> 


and  computes  the  Fourier  transform  of  the  LSF.  The 
diffraction  limited  Fourier  transform  of  the  optics  is  now 
computed  if  the  MIF  of  the  optics  has  not  been  previously 
measured  experimentally.  Next,  the  program  finds  the  MTF  of 
the  atmosphere  by  dividing  the  Fourier  transform  of  the  LSF 
by  the  Fourier  transform  of  the  optics.  Finally,  by  curve 
fitting,  the  program  computes  a  single  value  for  Cn2. 

The  program  now  predicts  the  size  of  a  laser  spot  on  a 
target  using  the  calculated  value  of  Cn2.  This  value  of  Cn2 
is  used  to  calculate  an  MTF  of  the  atmosphere  which  is  then 
multiplied  by  the  Fourier  transform  of  a  source,  and  the 
Fourier  transform  of  the  optics.  The  program  then 
calculates  the  inverse  Fourier  transform  of  the  products  cf 
the  above  and  uses  the  Abel  transform  to  give  the  angular 
point  spread  intensity  distrioution.  From  this  data,  the 
fraction  of  energy  as  a  function  of  the  total  energy  within 
a  given  radius  R  is  calculated. 


II 


•  THEORETICAL  Z3NSI3ERATI0NS 

The  theory  of  laser  beam  propagation  through  a  -urbulent 
medium  has  been  explained  by  Crittenden,  and  others,  and  is 
re-emphasized  here  for  continuity  purposes  [Ref.  2].  Since 
the  effects  of  turbulence  on  laser  weapons  is  of  major 
concern,  measurement  and  prediction  capability  for  these 
effects  or.  laser  beams  is  entirely  relevant. 

Due  to  the  existence  of  reciprocity,  these  ideas  apply 
to  either  laser  designators  or  imagers  [Ref.  5].  In  the 
study  of  the  theoretical  model  by  Fried,  the  effects  of 
atmospheric  turbulence  are  investigated  [Ref.  6  ].  This 
model  uses  the  idea  of  a  long  term  optical  transfer  function 
(OTF)  when  considering  atmospheric  turbulence.  The  long 
term  OT?  results  from  simply  taring  an  image  of  sufficiently 
long  term  which  sees  effectively  all  possible  turbulence 
configurations. 

When  observing  tae  effect  of  diffraction  as  the  result 
of  the  finite  aperture  of  a  point  source,  it  can  be  seen 
that  the  resulting  image  is  not  uniform.  Considering  the 
point  spread  function  to  have  the  same  shape  in  the  image 
plane  regardless  of  its  position,  the  image  function  is  the 
result  of  the  convolution  of  the  source  function  and  the 


op-tics  diffraction  function.  The  convolution  theorem,  as 
described  in  Fourisr  transfoci  theory,  yields 

i(v(x),v(y))«H  (v(x),v  (y))  *o(v(x)  ,v(y))*M(v(x),v(y))  (2.1) 

where 

i(v(x),v(y))  *  Fourier  transfoci  of  the  image  function 

H  (v  (x)  ,v  (y)  )=  Fourier  transform  of  the  ootics  diffraction 
function 

o(v(x),v(y))  =  Fourisr  transfoci  of  the  object  function 

M  (v  (x)  ,  ( v ( y ) )  *  Modulation  Transfer  Function  of 

the  atmosphers 

v  *  spatial  frequency 

A  point  source  such  as  a  Laser  can  be  analyzed  in  two 
dimensions  using  a  point  spread  function.  Through  Fourier 
transform  theory,  the  image  point  spread  function  is 
transformed  into  a  two  dimensional  OTF  of  the  optical 
signal.  This  probLem  may  be  simplified  by  scanning  the 
image  point  spread  function  using  a  vidicon  or  slit-scanning 
system.  In  applying  the  convolution  theorem,  the  Fourier 
transform  of  the  image  point  spread  function  is  multiplied 
with  the  Fourier  transform  of  the  optical  system  resulting 
in  the  Fourier  transform  of  the  overall  system.  The  Abel 


transfers  describe!  by  Griem,  is  appliei  to  this  rasa!-,  to 
re-transform  the  oie-dimensional  image  LSF  into  a  two 
dimensional  image  PSP  [Ref.  7]. 

As  demonstrated  by  Crittenden,  and  others,  a  numerical 
value  for  C„*  may  be  obtained  by  curve  fitting  using  the 
following: 

M  *  exp(-21.49*Cn**Z*f5/3*A"1/3)  (2.2) 

where 

8  *  MTF  of  the  atmosphere 
Z  *  Range  in  meters 

f  *  f*v  »  Angular  spatial  freguency  in  cyclss/radian 

F  =  focal  length  of  the  optical  system 

v  a  linear  spatial  frequency  in  cycles/meter 

A  *  wavelength  in  meters 

Cn2  is  obtained  by  a  linear  regression  of 

ln(B)  »  -21.49<Cn2*Z*fs/3*A“l /3  (2.3) 

where  Cn*  is  the  only  parameter. 
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III.  EXPERIMENTAL  procedure 


The  experiment  is  performs!  using  two  different  laser 
sources,  helium-neon  (He-Ns)  and  gallium-arsenide  (Sa-As). 
Figure  3.1  shows  a  block  diagram  of  the  experimental  sat-up 
used  for  both  lasers. 

The  measurement  of  Cn*  along  the  optical  path  is  made  by 
using  a  via  icon  and  telescope  a;  the  far  sad  cf  a  he  corridor 
in  the  basement  of  Spanagal  Hall.  Atmospheric  turbulence  is 
produced  by  nine  overhead  hot  air  ventilators.  The  optical 
equipment  similar  to  that  described  in  refarence  3,  consists 
of  a  6  inch  diameter  Cassegrain  telescope  with  a  90  inch 
focal  length  [  Ref  •  3]. 

Using  neutral  density  absorption  filters  to  attenuate 
the  intensity  of  the  laser  beam,  the  telescope  is 
illuminated  at  its  input  aperture.  The  light  image  from  the 
telescope  is  split  by  a  beam  splitter  with  one  beam  sent  to 
the  vidicon  and  the  other  to  a  slit  scanner  for  comparison 
of  the  two  systems.  The  light  is  then  transformed  to  an 
analog  signal  and  recorded  on  a  Panasonic  (NV-9240)  video 
tape  recorder  for  processing.  The  lata  from  the  slit 
scanning  system  are  recorded  on  a  precision  instrument  tape 
recorder  for  analysis.  After  recording  the  video  signa]. 
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analysis  is  accomplished  using  a  Tektronix  468  Digital 
Storage  Oscilloscope  and  a  Hewlett  Packard  9825  compu ting 
System.  The  linearity  of  the  Panasonic  tape  recorder  was 
demonstrated  in  Crager's  thesis.  This  is  also  confirmed  by 
observing  the  images  as  real  time  and  recorded 
displays  [Ref.  3]. 

The  468  oscilloscope  is  set-up  in  the  following  manner. 
An  input  signal  to  channel  A  is  located  on  the  scope  using 
the  non-storage  mode.  The  signal  is  then  displayed  using 
the  A  INTENSITY  switch  of  the  horizontal  display.  This 
intensified  zone  is  used  to  position  the  3  sweep  (delayed) 
to  the  desired  location  within  the  A  sweep  interval  to 
obtain  an  expanded  view  of  a  waveform  for  examination.  Once 
the  waveform  is  centered  on  the  scope,  the  horizontal 
display  is  switched  to  B  DLY'D.  This  is  done  to  facilitate 
the  digital  storage  circuitry  time  base  by  using  the  setting 
of  the  B  TIME/DIV  switch. 

The  waveform  is  expanded  in  time  by  decreasing  the  A 
TIME/DIV  switch  setting  and  moving  the  waveform  back  to  the 
middle  of  the  scope  using  the  delay  time  position  control 
dial.  When  a  representative  waveform  is  obtained,  the  B 
TIM  E/D  17  switch  is  then  used  in  conjunction  with  the  delay 
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time  position  control  dial  to  expand  the  tins  scale  until 
the  single  central  maximum  of  the  waveform  is  centered  and 
one  horizontal  sweep  is  displayed. 

Once  this  waveform  is  satisfactorily  obtained,  the  AVG 
storage  mode  is  selected.  The  468  will  average  the  input 
signal  for  a  selected  number  of  sweeps  and  display  the 
accumulated  waveform.  All  data  for  this  thesis  use  32 
sweeps  for  ea^h  average.  The  468  is  now  ready  to  transfer 
data  when  interrogated  by  the  HP  9825. 

The  computer  is  the  controller  for  all  interfacing 
operations  with  only  a  minimal  amount  of  operator 
interaction.  The  operator  interface  is  mainly  to  ensure 
that  the  equipment  is  properly  set  up  and  to  select  if  data 
are  to  be  plotted.  Digitization  of  the  analog  signal  is 
accomplished  by  the  Tektronix  45^.  Processed  waveform  data 
are  transferred  from  the  microprocessor  memory  to  the 
Storage  Display  RAH  [Ref.  8]. 

The  controlling  program  of  the  HP  9825  interrogates  the 
468  via  the  IEEE  488  interface  bus.  The  468  receives  the 
data  request  from  the  controller  and  sends  the  waveform 
message,  both  preamble  and  data.  The  waveform  message  is 


stored  in  the  calculator  memory  for  further  processing 


When  the  message  is  completed,  the  468  concludes  with  an  end 
of  instruction  recall ator  and  the  controller  takes  ccrrrol 
of  the  bus  again. 

Data  processing  of  the  digitized  signal  begins  with  the 
data  being  stored  on  magnetic  tape  for  further  use  by  the 
computer.  The  data  are  stored,  processed  and  plotted  by  the 
main  program.  Subroutines  are  called  as  necessary  for  their 
specific  uses.  The  two  signal  waveforms  are  recorded  and 
processed.  First,  the  signal  from  a  laser  beam  incident  at 
the  aperture  of  the  telescope  is  recorded.  Figures  3.2  and 
3.3  show  these  recorded  data  wiich  are  referred  to  as  the 
source  data.  .Second,  for  calibration,  the  signal  from  a 
laser  beam  with  a  diffraction  grating  in  place  at  the 
aperture  of  the  telescope  is  recorded.  The  entire  system  is 
calibrated  by  using  a  grating  that  consists  of  closely 
spaced  vertical  bars  in  front  of  the  telescope.  This 
produces  a  diffraction  pattern  in  the  image  plane  and  allows 
for  the  calculation  of  the  scale  factor.  These  recorded 
data  are  referred  to  as  the  sca-le  data  and  are  plotted  in 
Figures  3.4  and  3.5. 

During  the  recording  process,  the  signal  from  the  laser 
beam  is  recorded  on  video  tape  while  the  signal  is  observed 


Source  Data  for  Ga-As  Laser 
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Fig.  3.5  Plot  of  Scale  Data  for  Ga-As  Laser 


on  the  45  3  oscilloscope.  This  procedure  is  important  in 
data  taking  since  laser  alignment  and  telescope  focusing  are 
often  very  sensitive  to  minor  movements  of  either.  Upon 
completion  of  recording  the  data  on  video  tape,  the 
experiment  is  concerned  with  the  compilation  and  analysis  of 
data.  The  first  item  to  be  determined  is  the  scale  factor. 
Once  the  scale  factor  has  bean  determined,  it  will  remain 
constant  throughout  the  calculations,  unless  the  telescope 
focal  length  is  changed. 

The  scale  factor  is  calculated  in  the  following  manner. 
The  number  of  points  between  the  peaks  in  the  plot  of  the 
scale  data  (diffraction  grating  in  place  in  front  of  the 
telescope  aperture)  is  measured.  When  the  distance  bet  wee- 
peaks,  the  spacing  between  the  bars  in  the  grating,  and  the 
wavelength  of  the  laser  are  known,  the  scale  factor  can  be 
calculated  from  the  relation 


sin  (theta)  =  X /d 


(3. 1) 


sin  (theta)  ~  theta 


(small  angle  approximation) 


(3.2) 


sf  *  (theta)  /ndp 


(3.3) 


sf  »  A / (d*ndo) 


(3.4) 


where 


sf  =  scale  factor  in  radians  par  point 
A  =  wavelength  in  meters 

d  =  spacing  between  lines  on  grating  in  asters 

ndp  =  numbep  of  points  between  rsntral  maximum 
and  first  order  diffraction  peak 

This  standardizes  the  data  for  the  abscissa  of  the  plots  in 
rad-ians  per  point. 

The  program  now  takes  the  transferred  data  and  computes 
a  point  spread  function  as  shown  in  Figures  3.5  and  3.7. 
Figures  3.8  and  3.9  represent  the  point  spread  function 
after  integration  using  equation  (1.1)  to  obtain  a  line 
spread  function.  Next,  the  Fourier  transform  of  the  line 
spread  function  is  calculated.  These  curves  are  plotted  in 
Figures  3.10  and  3.11. 

The  diffraction  of  the  optics  is  computed  and  plotted  as 
shown  in  Figures  3.12  and  3.13.  Following  the  same  method 
as  before,  the  line  spread  function  of  the  optics  function 
is  calculated  and  pLotted  in  Figures  3.14  and  3.15.  The 
program  now  takes  the  Fourier  transform  of  the  optics 


function  LSF 


These  data  are  plotted  in  Figures  3.16  and 


Fig.  3.8  Plot  of  Point  Spread  Function  for  He-Ne  Laser 


Fig.  3.7  Plot  of  Point  Spread  Function  for  Ga-As  Laser 


3.  17 


The  Fourisr  transforms  of 


the  system  ana  cp 


lie 


divided  point  by  point  to  yiaii  the  MTF  of  the  atmosphere, 
as  shown  in  Figures  3.18  and  3.19.  The  MTF  of  the 
atmosphere  will  be  used  to  calculate  the  long-term  value  of 

C„2* 

Cn2  is  calculated  by  linear  regression  of  In  (MTF)  versus 
f s/3.  The  slope  of  the  regression  is  proportional  to  Cn2. 
Twc  simultaneous  equations  are  solved  using  Cramer's  rule 
and  the  matrix  ROM  of  the  HP  9825  .  The  equations  used  are 

a  *X  +  b*I  =  Y  (3.5) 

a  *X2  ♦  b*X  =  Y*X  (3.6) 

where 

a  =  slope  of  the  curve 
b  =  intercept 

X  =  summation  of  (I*sf )  **( 5/3) 

X**2  =  summation  of  (I  *sf)  **  (1 0/3) 

I  =  point  number  (total  number  =  256) 

Y  =  summation  of  natural  logarithm  of  MTF 
Y*X  *  summation  of  products  of  two  values 


Cn2  is  obtained  from  the  above  Information  using  equation 
(2.3).  This  will  yield 

Cn2  =  a/(-21.  49*Z*f**  (5/3)  *  **(-1/3))  (3.7) 

The  program  now  goes  to  the  prediction  phase  after 
calculating  C  n  2  foe  the  atmosphere.  If  a  Gaussian 
distribution  for  an  input  source  and  a  value  of  the  standard 
deviation  are  assumed  or  known,  the  resulting  source 
function  can  be  calculated.  A  plot  of  the  computed  source 
is  shown  in  Figure  3.20.  It  can  be  seen  from  this  plot  that 
the  half-width  at  half-maximum  is  approximately  4 
microradians.  As  before,  thja  line  spread  function  of  the 
computed  source  is  calculated  and  plotted.  This  is  shown  ir. 
Figure  3.21.  The  Fourier  transform  of  this  data  is 
calculated  and  plotted  in  Figure  3.22. 

The  program  now  multiplies  the  Fourier  transform  of  the 
computed  source  with  the  transform  of  the  system  (including 
the  atmosphere)  and  plots  the  result  in  Figures  3.23  and 
3.24.  Plots  of  the  inverse  Fourier  transform  are  shown  in 
Figures  3.25  and  3.25.  Next,  the  Abel  transform  is  computed 
and  plotted  in  Figures  3.27  and  3.23.  The  Abel  transform, 
as  described  in  Chapter  II,  transforms  a  one-dimensional 
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line  spread  function  into  a  two-dimensional  point  spread 
function.  Finally,  tae  fraction  of  power  inside  a  circle  of 
radius  R  is  calculated  and  plotted  as  shown  in  Figures  3.29 
and  3.30.  This  is  the  fraction  of  power  that  one  would 
expect  to  be  incident  on  a  target  using  the  measured  value 
of  atmospheric  turbulence  as  an  input  parameter.  A  detailed 
analysis  of  the  computer  program  is  the  subject  of  the  next 


3.8  Plot  of  Line  Spread  Function  for  He-Ne  Laser 


Fig.  3.9  Plot  of  Line  Spread  Function  for  Go-As  Laser 


3. 10  Plot  of  Fourier  Transform  of  LSF  for  He-Ne  Laser 


Fig.  3. 12  Plot  of  Optics  Diffraction  for  Hfc-Na  Laser 


Fig.  3. 13  Plot  of  Optics  Diffraction  for  Gcj-A9  Laser 


Fig.  3. 14  Plot  of  LSF  of  Optics  Function  for  He -No  Loser 


cs  Function  for  Go- As  Laser 


Fig.  3. 18  Plot  of  Fourier  Transform  of  LSF  of  Optics  Function  for  He-Ne  Laser 


Fig.  3.  17  Plot  of  Fourier  Transform  of  LSF  of  Optics  Function  for  Ga-As  Laser 


Atmosphere  for  He-Ne  Laser 
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Fig.  3. 19  Plot  of  MTF  of  Atmosphere  for  Ga-As  Laser 


3.20  Plot  of  Computed  Source 
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Fig.  3.22  Plot  of  Fourier  Transform  of  LSF  of  Computed  Source 


Fig.  3.25  Plot  of  Inverse  Fourier  Transform  for  He-Ne  Laser 


Fig.  3. 26  Plot  of  Inverse  Fourier  Transform  for  Ga-As  Laser 
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Fig.  3.28  Plot  of  Abel  Transform  for  Ga-As  Laser 


3.29  Plot  of  Power  Fraction  Inside  Circle  of  Radius  R  for  He-Ne  Laser 


Fig.  3.30  Plot  of  Power  Fraction  Inside  Circle  of  Radius  R  for  Ga-As  Laser 


IV.  COMPUTER  PROGRAM  IS  ALGORITHMIC  LARGHAGd 


The  program  presented  here  is  written  in  rhe  algorithmic 
language  as  outlined  by  Graham  [Ref.  9].  It  is  intended 
that  the  language  in  use  here  be  machine  independent. 
Comments  appear  in  parentheses. 


Alger  i  th  L  aser 
(Input  which  laser  is  used.) 

Input  Laser_Wavelength 
ESTRING  <«=  Type_o£_Laser 

(Input  if  the  particular  laser  source  data  are  recorded. 
The  following  lines  determine  if  the  laser  source  data  are 
already  recorded.  If  the  source  data  are  not  recorded  rhsn 
the  program  goes  to  the  transfer  subroutine  to  allow 
interfacing  with  the  468  to  effect  a  data  transfer.  The 
parameters  passed  to  the  Plot  subroutine  have  the  following 
significance:  the  first  parameter  determines  the  horizontal 

index,  for  example  512  points  or  the  number  of  microradians. 
The  second  parameter  tells  the  computer  which  file  the  data 
are  recorded  on  so  the  data  can  be  loaded  iQto  memory.  The 
final  parameter  determines  the  labeling  of  the  horizontal 


or 


a xis  e irner  ir.  nata  points,  arcrcradians, 


nr.es  per 


microradian .) 

If  Laser,Source_no t_recorded  then 
Call  Transfer  (Source, File) 

Call  Plot  (512  , Source, File,  0) 

End  If 

(Input  if  the  particular  laser  scale  data  are  recorded. 
Scale  data  are  recorded  with  the  diffraction  grating  in 
place.  If  the  scale  data  need  to  be  recorded,  the 


subroutine  Transfer  is  called  and  the  scale  data  are  stored 
on  a  separate  file.) 

If  Laser  Scale  not  recorded  then 


Call  Transfer  (Scale_File) 


Call  Plot  (512  ,Scale_fiLe,0) 


End  If 


Input  Scale_Factor 


(Calculation  of  the  scale  factor  is  detailed  in  Chapter  III. 
Comparison  of  the  time  base  of  the  recorded  data  is 
necessary  to  ensure  that  the  data  recorded  have  the  same 


time  scale  on  the  oscilloscope.  If  the  time  scales  are 
different  the  program  stops  and  displays  an  error  message. 
k  file  is  loaded,  the  waveform  preamble  is  searched  for  the 


"X  UNITS"  or  "XINCR"  oscilloscope  settings  and  -he 
particular  data  are  stored.  This  is  done  for  each  of  the  3 
recorded  files:  source,  scale,  and  collimated  data. 
Collimated  data  refer  to  the  same  laser  recorded  under 
controlled  conditions.  The  laser  beam  is  projected  in  a 
collimator  and  recorded  on  tape.  This  data  could  be  used  to 
remove  any  non-linearities  of  the  video  recorder,  electronic 
equipment,  or  receiving  optics.  Data  elements  are  compared 
to  ensure  consistency  of  inputs.) 

Input  Laser_sourre 

B<==  po  s(Laser_Sou  rce,  » R:  • ) 

R1<==  v  al  ( *  B+2*  ,3*  4') 

3<  ==  po  s(Laser_sourcs, 'IT: ) 

ASTRING<=  =  Laser,  3  ource  ('  B+3  '  ,  '3+4') 

Input  ( Laser  _Scale< 

B<==  po  s  (Laser,  Sc  a  le, '  R : ' ) 

R2<==  val  { '  B+2'  ,  3*  4 ' ) 

B<==  po  s(Laser_Scale, '  IT: ') 

BSTRING<==  Laser  ,3  cale.(  'B*3' ,  'B+4») 

Input (Colli mated, 3  ource) 

B<=*  po s (Collimate d_So  urce, '  R:  ') 

R3<=*  v al ( ' B+2' , ' 3  +4) 

B<®=  pos(Collimate  d_So  urce, '  IT:  ') 
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CS?RING<*=  Collimated^ Source  (,5  +  3,f’3  +  ti') 

If  R  1=R  2=R 3  then. 

If  ASTRING  =  BSTRING  =  CSTRING  then 
Input  RD ATA (File.* umber) 

End  If 

Else 

Output  * ERR3R  IN  TIME  REALS' 

Stop 
End  If 

A<==B<  =  =C<=*=0 

(Array  is  initialized  to  0) 

Do  for  I  =  1  to  5 1 2 

YD  AT  A  (I)  <=-  R  DATA  (I) 

End  Do 

(This  section  takes  the  data  from  the  waveform  message  as 
transferred  from  tha  468  and  removes  any  DC  background  from 
an  image  as  well  as  subtracting  the  minimum  value  in  order 
to  "zero"  the  data.  This  makes  y=Q  as  the  starting  point 
for  both  calculation  and  plotting  purposes.  Data  up  to  51 
and  greater  than  463  are  summed  and  stored  in  B. 

This  represents  the  "wings"  of  the  curve.  This 
process  removes  the  DC  background  from  an  image.) 

Do  for  I  *  1  to  51 2 
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A<=*  A+YDATA  ( I) 


If  I  <  51  then 


B<«  YDATA(I)  ♦  B 


End  If 


If  I  >  463  then 


B<==  YDATA(I)  ♦  B 


End  If 


End  Do 


(Output  for  debugging) 


Output  A,  B 

(Average  of  the  sum  of  data  values  in  the  "wings".) 

D<==  B/100 

(This  is  the  total  sum  of  all  data  less  ths  "wings"  value.) 
k<-~  A-512*D 

(All  data  are  summed  and  stored  in  C.  Each  data 
is  compared  with  the  previous  total  divided  by  2. 

This  will  find  the  index  value,  I,  for  which 
the  maximum  value  exists.) 

Do  for  I  *  1  to  512 

YDATA(I)<  =  =  Y  DATA  ( I)  -  D 
C<»»  YDATA  (I|  ♦  E 


If  C  <  A/2  then 


J<=*  I 


p.d  If 


End  Do 

(Output  for  debugging.) 

Output  C,J 

(The  curve  is  shifted  so  that  the  maximum  value  starts 
at  the  origin  of  the  horizontal  axis.  This  is  called 
the  point  spread  funotion  and  is  recorded  as  RDATA.) 
Do  for  K  =  1  to  256 
I  DAT  A  (K)  <=  = 

YDATA  (  (J-K  +  1)  mod51 2  +  1 1  ♦  YDATA  ( ( K**  J  —  1)  mod512«-1) 

End  Do 

Do  for  I  =  257  to  512 
RDATA  (I)  <  =  =  0 
End  Do 

Do  for  I  =  1  to  256 

RDATA  (!)<==*  I  DATA  (I) 

End  Do 

(Determines  if  tabular  or  plotted  data  are  desired.) 

If  PSF _ OUT P Unwanted  then 
Call  Table  (RDATA) 

End  If 

If  PSF_PLOT_wanted  than 

Call  Plot  (32,  RDATA,  1) 


[♦] 


(Calculates  the  lias  spread  function  from  the  point 


spread  function.) 

Call  LSF 

(Determines  if  tabular  or  plotted  data  are  desired.) 
If  LSF_OUTPUT_wanted  then 
Call  Table  (RD  AT  A) 

Else  If  LSF_PLOT_v anted  than 
Call  Plot  (32,  RDATA,  1) 

End  If 

(Calculation  of  Fourier  Transform  of  LSF.) 

Input  RDATA 

(The  maximum  value  of  the  array  IDATA  is  at  1*1.  This 
is  put  back  into  array  RDATA  and  folded  over  so  that 
IDATA  (257)  to  IDATA(5  12)  is  the  mirror  image  of 
IDATA  (1)  to  IDATA (255).  Then  this  array  is  recorded 
as  RDATA  and  the  Fourier  transform  calculated.) 

Do  for  I  *  257  to  512 
RDATA  (I)<*=  0 
End  Do 

Do  for  I  *  1  to  256 

RDATA  (514-1)  <«  RDATA  (I) 


(Determines  if  tabular  or  plotted  data  are  desired.) 

If  FXFORM_OF_L5F_3  DTP3T_if  ant  ed  then 
Call  Table  (RD  AT A| 

Else  If  ?XFORti_OF_LSF_PLOr_wanted  then 
Call  Plot  (32,  RDATA,2> 

End  If 

(Dia  of  Obscur/Dia  of  Obj  Lens  in  meters, RATIO*. 064/. 164) 
(This  is  the  measurement  of  the  ratio  of 
diameters  of  obscuration  to  telescope.) 

(Hi*  Scale  of  data  in  microradians=1. 60) 

(Diameter  of  telescope  Obj  Lens  in  meters , OBJECT* 

0.  164  meters) 

If  Laser_Sourca  =  He-Se  then 
Wa velength<==  6.328*a-7 

Else 

Wa velength<==  9.05*e-7 
End  If 

(This  begins  the  calculation  of  the  Fourier  transform 
of  the  optics  using  the  Airy  function.  All  of  the 
data  regarding  the  optics  are  recorded  on  a  separate  file.) 
(These  next  four  statements  fir  the  constants  in  the 


argument  of  the  optics  function.) 


R1<==  Scale  Factor 


D<==  RATI0**2 


H<==  1-D 


Z<==  PI  *Scale_Fact  or  *D  bject  *  1 *e-6/Wavelength 


Do  for  I  *  1  to  256 


Y<  =  Z*(I-1| 


If  Y  >  30  than 


RDATA  (I|  <==  0 


Else 


RDAIA (I)  <==  (AIRY (Y) -D*AIRY (Y*R  ATIO) ) **2/3 **2 


End  If 


K<=  I 


End  Do 


(Record  RDATA  on  separate  fiLe) 


(Determines  if  a  plot  of  the  optics  function  is  desired.) 
If  0PTICS_F0NCTI3N_PL3T  wanted  then 
Call  Plot  (32,  RDATA,  1) 

End  If 

(Calculates  the  line  spread  function  of  the  optics 
function. ) 

Call  LSF 

(Determines  if  a  plot  of  the  LSF  of  the  optics  function  is 


If  PLOT _OF_LSF_OPTICS  wants}  then 
Call  Plot  (32,  RDAT  A,  1) 

End  If 

Do  for  I  =  257  to  512 
RDATA  (I)  <==  3 
End  Do 

it  a  are  folded  over  again  ana  tie  Fourier 

calculatsd .) 

Do  for  I  =  1  to  256 


RDATA  (514-1)  <==  RDATA(I) 

End  Do 

RDATA  (2  57)  <-*  BOAT  A  (256) 

(Record  RDATA) 

(Fourier  transform  of  the  LSF  of  Optics  is  calculatsd.) 
Call  FXFORM 

(Detsrmines  if  plot  of  Fourisr  transform  of  Optics 
function  is  desirsd.) 

If  PL0T_0F_FXF0R1_0F_3PTICS  wanted  then 
Call  Plot  (255  , RDATA,  2) 

End  If 

(Calculation  of  tha  quotient  of  two  Fourier  transforms.) 
(This  section  takes  the  Fourisr  transform  of  the  system 
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and  divides  it,  point  by  point,  by  the  Fourier  transform 
of  the  optics.  This  yields  tie  modulation  transfer 
function  (MTF)  of  the  atmospiere.) 

Input  SDAr&_OF-FXFORM_OF_SO[JEC2 
Do  for  I  =  257  to  512 
RDATA  (I)  <  =  =  3 
End  Do 

Do  for  I  =  1  to  256 

RDATA  (514-1)  <==  RDAIA  (I) 

End  Do 

Input  R  DAI A_OF_FX  F  ORM^OF^DPI ICS 
Do  for  I  =  1  to  256 

RDATA(I)<  =  =  IDATA  (I) /RDATA  (I) 

End  Do 

(Determines  if  a  plot  of  the  3TF  is  desired.) 

If  PL3T_0F_MTF_0F _ ATM3S_ wanted  then 
Call  Plot  (32,  RDAIA, 1) 

End  If 

(Calculation  of  Cn2  by  curve  fitting. 

This  procedure  is  described  in  Chapter  III.) 

(All  "8"  variables  used  in  this  calculation  are 
initialized  to  0.) 


R17<»=a  18<  =  =  R1 9<  =  = R2Q<=  =  R21<  =  =  R22<==R23<=  =  R24<==R25<==0 


(Analysis  of  data  shows  that  computer  tins  is  wasted 
beyond  96  points.) 

Do  for  I  =  1  to  95 

(R18  -  when  summed  is  the  (f(i»)l°/3.  This  is  the 
total  of  the  angular  spatial  frequencies  squared.) 

R19<==  ((I-1l  *Scale_Faotor)  *♦(10/3) 

R  1  70  =  R17  «■  R18 

(R19  -  when  summed  is  the  (£(ii)5/3.  This  is 
the  total  of  the  angular  spatial  f requenoies. ) 

R  1  90=  ((1-1)  *Scale_Faotor)  **  (5/3) 

R20O=  HI 9  R20 

(R21  -  is  the  total  number  of  points.) 

R2  10=  I 

Do  While  RDAT  A  (It  >  3 

(R22  -  when  summed  is  the  total  of  the  product  of  the 
natural  logarithm  of  each  point  with  R19.) 

S22<==  1  n  (RD  ATA  (II  )  *R19 
?,  230=  R  22  R23 

(R2h  -  when  summed  is  the  total  of  the  natural  logarithm 
of  each  point.) 

R  240=  1  n  (RD  ATA  (I)  ) 

R25<==  R  24  ♦  R25 
End  Do 


5  5 


DO 


End 

(Output  matrix  values  for  deougging.*- 
Output  R17, R20, H21 ,R23, R25 
(Assign  value  to  matrix. 

This  section  takes  the  above  oalculated  values  and 
sets  up  a  matrix.  The  matrix  equation  solved  is 
B  =  A-‘*C.) 

(Create  3  matrices  A  (2,  2)  ,3(2,1)  ,C(2,1).) 

A  (1,1)  <=*  R20 
A  ( 1  r  2)  <*=  R 21 
A  (2,  1)  <==  R 17 
A  (2,  2)  <=*  R  20 

0 

INVERT^  MAT  RIX_A 
C  (1,  1)  <—  R 25 
C  (2,  1)  <  — R2 3 
B<  ==  A  *  C 

(Output  matrix  values  for  debugging.) 

Output  B  (1  ,  1)  ,B  (1,  2) 

(Calculation  of  Cn2  using  equation  (2.  3) . ) 

* 

R22<==  B (1 »  1)/(-21 . 49* Ranged R 20*  (Wavelength**  (-.  3333  3) ) ) 
Output  ,CNSQ*,fR22 

(The  program  now  starts  the  prediction  phase.  It 
calculates  a  source  assumed  to  have  a  Saussian 
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distribution  by  mains  of  tha  computed  source  pattern 
below.  Using  similar  Fourier  transform  techniques 
the  program  uses  the  calculated  value  of  Cn2  and  predicts 
the  power  incident  on  the  target.) 

(Computed  Source  Pattecn:A=AD*exp  (-x**2/2*Sigma**2) 

(This  is  an  "arbitrary"  Gaussian  source  pattern 
with  a  standard  deviation  for  Sigma  set  equal  tc  2.) 

AC<==  A  _3  2  R  0 

C<==  S igma*Scale_F actor 

Do  for  I  *  1  to  25  6 

F<=  (1-1)  *Soale_Factor 
3  <  =  F**2/  (2*C**2) 

Do  While  3  >  13 

HD  AT  A  (I)  <==  0 
End  Do 

RDATA(I)<=  =  A3*exp  (-G) 

End  Do 

(Record  BOAT A.  This  records  the  computed  source.) 
(Determines  if  plot  of  computed  source  is  desired.) 

If  PLOTJ3F_COMPUIED_SDURCE_*anted  then 
Call  Plot  (32,  R  DATA  ,  0) 

End  If 

(Calculates  the  line  spread  function  of  the  computed 
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source. ) 

Cali  LSF 

(Determines  if  a  plot  of  the  LSF  of  the  computed  source 
is  desired.) 

If  PLDT_0F_LSF_0F_C0MP_S0[JR:E  =  1  then 
Call  Plot  (32,  RDATA , 0) 

End  If 

(Calculation  of  the  Fourier  transform  of  LSF  of  computed 

source. ) 

Do  for  I  =  257  to  512 
RD  AT  A  (I)  <*  =  0 

End  Do 

Do  for  I  *  1  to  256 

RD  ATA  (514-1)  <==  RDAIA(I) 

End  Do 

RDATA(2  57<==  RD  AT  A  (256 ) 

Call  FX  FORM 

(Determines  if  a  plot  of  the  Fourier  transform  of  the 
LSF  of  the  computed  source  is  desired.) 

If  PLOT_OF_FXFDRH_OF JTOMParSD_SO(JRCE_wanted  then 
Call  Plot  (255  ,PDATA,2> 

End  If 

(Calculation  of  the  product  of  two  Fourier  transforms.) 
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(This  is  the  Fourier  transform  of  the  computed  scarce 
with  the  STF  cf  the  atmosphere.) 

Do  for  I  =  1  to  256 


RDATA  (I)  <=  = 

RD  AT  A  (I)_CO)iPUTED_SDtJRDE*RDATA  (I)  _FXFK_OF_^TF 
End  Do 


(Product  of  FXF0RS3  of  Source*Atmosphere*Optics.) 

RD AIA (!)<== 

RDATA( 1 ) _PRODUCT_D F_2_FXFMS“ RD ATA ( 1 )  _FXFM_OF_OPTICS 
Do  for  I  =  2  to  2  56 

RD  ATA  (I)  <=  = 

3D  ATA (I) _PRDD(JCT_0F_2_FXFKS*RDATA (I)  _FX F«_OF_DPTICS 
RDATA  (514-1)  <==RDATA(l! 

End  Do 


RDATA  (2  57)  <==  RDATA  (256) 

(The  above  data  predicts  what  the  Fourier  transform 
of  the  entire  system  is,  using  the  calculated 
value  of  Cn2  for  the  atmospheric  turbulence.) 

(Determines  if  a  plot  of  the  result  of  Fourier  transform 
is  desired.) 

If  PLOT_OF_FXFORH_  PRODUCTS _w  anted  then 
Da  11  Plot  (32,  RDATA, 2) 

End  If 
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(Inverse  Fourier  transform  gives  target  LSF.) 

Call  IN  VERSE_FXF0  8  M 

(Determines  if  a  plot  of  the  inverse  Fourier  transform 
is  desired.) 

If  PLOT_OF_INVERSE_FXFROM_vanted  then 
Call  Plot (32,  RDAT  A, 0| 

End  If 

(Converts  target*s  one  dimensional  LSF  to  a 
two  dimensional  PSF  by  Abel  transform.) 

Call  ABEL 

(Determines  if  a  plot  of  the  Abel  transform  is  desired.) 

If  PLOT_OF_ABEL_XFORM_wanted  then 
Call  Plot  (32,  RDAIA,0» 

End  If 

(Calculates  the  fraction  of  power  inside  circle  of  radius 
R.  This  predicts  the  fraction  of  power  that  will  be 
incident  on  target. ) 

RD  AT  A  ( 1  )  <**  0.25*PI*RDATA(1> 

Do  for  I  *  2  to  236 

RDATA  (!)<*=  2  *PI*  RD ATA  (I)  +RDATA(I-1) 

End  DO 

(Determines  if  a  plot  of  the  fraction  of  power  is  desired.) 
If  PLOT  OF  POWER  wanted  then 
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all  Plot  (32,  RDAr  A,  0) 


End  If 
End  Laser 

(The  following  subroutines  ire  used  in  the  Algorithm 
Laser.) 

Subroutine_Transfer( File. Number! 

(The  transfer  subroutine  gets  caw  data,  both  preamble 
and  data  from  the  waveform  message  sent  by  the 
468.  It  processes  the  message  by  finding  the  minimum 
value  of  the  array  then  subtracting  this  value  from 
each  element  in  the  array.  This  "zeros"  the  array.) 

Output  ‘Ensure  Equipment  set  up  properly' 

Output  'Continue  when  ready' 

3<==  (-1) 

Do  While  E  *  -1 

("DATA"  is  the  I/O  buffer  where  the  data  from  the  468  are 
sent.  The  status  of  the  buffer  is  read  while  data  are  being 
transferred.  When  the  transfer  is  complete,  the  interface 
is  cleared.) 

Data<«  Trans ferred.Data 
If  Data_Traas fer. Complete  then 
E<=*  0 
End  If 
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Clear  Interface 


N<==  1 

(This  is  where  healer  format  stops  and  lata  begin.) 

A<==  po  s(Data,  *  S' ) 

Do  for  I  =  A+2  to  687  by  16 
Do  for  J  =  1  to  1 6 

If  J*I>5  88  t  hen 
RJC  ==  3 
N<  =  =  N*1 
End  If 

R  J<ss  na  m(Data  (I  +  J)  ) 

(The  numerical  values  of  each  element  are  stored  in  YDATA.) 
YDATA  (S|  <==  RJ 
N<==  N*1 
End  Do 

End  Do 

E<==  minimum  (YD Ar A) 

Do  for  I  *  1  to  51  2 

YDATA  (I)  <=*  YDAI  A  (I)  -E 
RD  AT  A  (I)  <=*  I  DAT  A  (I) 

End  Do 

Store_RDATA  (Filers umber) 


Return 


End  Transfer 

(This  subroutine  calculates  ths  Line  Spread  Function 
for  the  previously  recorded  point  spread  function 
by  using  equation  (1.1). j 
Subrout ine_LSF 

(Depending  on  which  file  number  has  been  passed  down 
from  the  calling  subroutine,  either  scale  or  source 
data  are  loaded.) 

If  File_Number  =  Scale_Data  then 
Input  Scala_Data 
RD  ATA<=«SCALE_DAr  A 

Else 

Input  RDAI A 
End  If 

Do  for  I  *  1  to  51 2 
ID  AT  A  (I)  <==  D 
End  Do 

Do  for  I  *  1  to  256 

ID  ATA  (I)  <  =  =  3  DATA  (I) 

End  Do 

(Plots  have  shown  that  computer  time  is  lost  and  no 
valuable  information  is  gained  beyond  about  24  points.) 
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Do  for  I  =  1  to  24 


J<  —  1 
Output  I 
Q<=*  IDATA  (I| 

Do  While  R  <  24 

R<==  S2RT  (I*I  +  J*J) 

Q<== 

Q+2*(  (1-  fraction  (5)  )  *i  (ir.t  (3) )  ^fraction  (?)  *i  (int  (F)  *  1» 
J<==  jfl 
End  Do 

ID  AT  A  (I)  <==  2 
End  Do 

Do  for  I  =  1  to  2 '4 

RD  ATA  (I)  <==  I  DATA  (I) 

End  Do 

If  File_Number  =  Scale_Data  then 
Store_Scale_D  ata 

Else 

Stoce_RDATA 
End  If 
Return 


End  LSP 


(This  subroutine  calculates  the  Fourier  transform 
of  the  given  data  using  the  Cooley-Tukey  Algorithm. 
If  the  inverse  statement  is  true,  then  the  inverse 
Fourier  transform  is  calculated.) 

Subrout in e_FXFORM 
Set_Radian_Mode 

(2**9=512  which  ie  the  number  of  points.) 

N<==  9 

Do  for  I  «  1  to  512 
ID  AT  A  (I)  <==  D 
End  Do 

If  File_Sumber  *  Scale^Data  then 
Input  Scale_3ata 

Else 

Input  RDAIA 
End  If 

T<=»  PI/2**(N-1) 

Do  for  J  *  3  to  2**  (3-1)  -1 
(BI  is  a  bit  inversion  subroutine.) 

Call  BI(J,P,N-1) 

C<*«  co  e  (P*I) 

If  IHVER SE_FXFORM  =  True  then 
Flg7<=»  1 
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End  If 

P<==  sia  (P*D  *(1-2*Flg7) 

Do  far  I  =  2*R0*J+1  to  2*R0*J+R0 

HK  ==  RDATA(I) 

R2C==  RDATA  (I  +  RO) 

R3<==  IDATA(I) 

R4<==  IDATA  (I ♦  RO) 

RD A  TA (I) <=»  3 1 ♦R2*C  +  R4*P 

IDATA  (I)  <==  33+R4*C-R2*P 

RDATA  (I  +  RO)  <  =  =  R1-R2*C-R'4*P 

♦ 

IDATA  (I  +  P.O)  <"  R3-R4*C+R2*P 
End  Do 
End  Do 
Output  M 
End  Do 

(This  section  is  for  re-oriering  the  block.) 
Do  for  I  =  0  to  2**N-1 
Call  BI(I,JrH) 

Do  While  I-J  <  or  *  3 


Do  While  I  i  J 


?<*=  RDATA(I«-1)  /SQRT(2**N) 


Z<  =  =  IDATA  (If  1)  /SQRT(2**N) 

RDATA  (1  +  1)  <  =  =  R DATA  (J  +  1 1 
IDATA  (J  +  1)  <  =  =  I  DAT  A  (J*  1 ) 

RDATA  (J ♦1)  <s=  p 
IDATA  (J+  1)  <  =  =  Z 
End  Do 

RDATA  (If  1)  <==  RDATA  (1+1) /SQRT  (2**N) 

IDATA  (I*  1)  <*=  ID  AT  A  (I  +  1)/SQRT  <2**N) 

End  Do 

If  File_Nuaber  =  Scale_Data  then 
Stora_Soale_Data  (File^Number) 

Else 

Store_3D  ATA(File_'l  a  a  bar) 

End  If 

End  Do 

Set_Degrees_Mode 
Return 
End  FXFORM 

(This  subroutine  takes  the  binary  number  PI  containing 
P3  bits  and  inverts  it  end  for  and,  e.g.  310111  become 
111010.) 

Subroutine  BI(P1,P2,P3) 

P2<3=  0 
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?4<=  =  PI 

Do  for  Z  =  1  to  ?3 
?  4  <==  P4/2 
?  2  <==  2*P2 

If  fraction(P4)  #  0  than 
P  2<==  P  2  ♦  1 
End  If 

P  4  <==  int(?'4| 

End  Do 
Return 
End  31 

(This  function  calculates  Airy ( x>  =2*J  1  (x|  /x  where  J1  (x) 
is  the  Bessel  function  cf  order  one.) 

Function  AIRY<P1| 

If  PI  <  0  then 

Output  ' ERROR -AR3  UMENI  LESS  THAN  3* 

Stop 
End  If 

If  PI  =  0  then 
R  4  <==  1 
Return  R4 
End  If 
R5<»a  0 
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If  PI  >  15  then 


?. 6  <=  =  90  +  PI  /2 
R  1  2<==  1. 4  *P 1  +  60/PI 
End  If 

If  PI  <  5  then 

R 6  <==  20  *■  1  3  *P  1- P  1**2/3 
R 1  2<==  6  «■  Pi 

Else 

R6<==  20  ♦  13 *Pl-?1**2/3 
R12<==  1. 4*P1  +  50/PI 
End  If 

R12<=*  maximum  (512>  ,int  (3+P1/4) ) 
Do  for  a  »  312  -.3  R6  by  3 
P  8<==  1*9-29 
R13<=  =  R14<==  0 
If  M/2  *  int  ( M/2)  than 
Flg13<==  0 

Else 

Flg10<==  1 
End  If 

Do  for  J  =  1  to  H-2 

R  1 5<==  2*(M-J)  *R9/P1-R13 
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R13<*»  R  8 


If  81—  J—  2  =  0  then 
R'4<==  R  15 
End  If 

If  FLglO  =  3  then 
Fig  10<==  1 

Else 

Fig 1Q<=  =  0 
Er.d  If 

R  14<  =  =  R  14+2*R8*Flg  10 
En  d  Do 

R1S<  =  =  2*R8/P  1-H1  3 
R14<  =  =  S14+R15 
H4  <==R4/R1  4 

If  (abs(P.4-R5)-abs<R’4*1*e-6))  <  or 
R4<==  2  *  R4/P 1 
End  If 
Return  R4 
R  5  <==  R4 
End  Do 

Output  'ACCURACY  SOT  OBTAINED' 

Return  R4 
End  AIRY 
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(The  subroutine  Abel  takas  a  oa a-dimensional 
line  spread  function  and  calcuL at en  a  two 
dimensional  point  spread  function.) 

Subroutine  ABEL  (Pila. Number) 

N<==  RDATA  (1) 

RDATA  (  1 )  <==  1.4+RD ATA( 1)  -1. 3  ♦‘RDATA  (2)  ♦.  4* RDATA  (3) 
Do  for  I  =  2  to  54 
M <  —  RDATA  (I| 

RDATA  (I)  <==  .  4*N+.2*!i-.  5  *RDATA  (1+1) 

N<=  M 
End  Do 

Do  for  I  =  1  to  54 

RD  ATA  (I)  <==  R  DATA  (I)  /(2*SQRT  (  (I  +  .  1)  * *2-1*1)  ) 
Do  for  J  =  1  to  54 
RDATA  (I|  <== 

RDATA  (I)  +RDATA  (J)  /SQRT  (  (J*.  1)  **2-1*1) 

End  Do 

RD  ATA  (I)  <*=  RDATA  (I) /PI 
Output  I 
End  Do 
Return 


End  ABEL 


(This  subroutine  prints  out  data  in  tabular  fori.) 
Subroutine_TABLE(Fil3_Nuaber) 

Input  RDAIA  (Filers  umber) 

Do  for  I  =  1  to  32 

Do  for  J  =  1  to  15 
Output  RDATA  (  16  (1-1)  ♦  J» 

End  Do 

Output  RDATA  (16*1) 

Return 
End  TABLE 

(This  subroutine  plots  the  desired  data.) 

(Each  time  a  plot  is  called,  tais  subroutine 
plots  the  particular  graph  and  labels  it. 

The  horizontal  axis  is  labeled  according  to  the  first 
and  last  parameters  passed  by  the  calling  subroutine.) 
Subroutine_PLOT  (Pi  ,Fi le_N umber ,  P3) 

Input  RDATA  (File_H  umber) 

A<  ==  0 
B<==  PI 

C<=»  minimum(RDArA) 

D<==  maximum  (RDATA ) 

XilIN<==  -.  1*  (B-A) 


XMAX<=  =  3+.05*(3-4) 


YMIN<=  =  C-.  1*  (D-C| 

YMAX<==  D+.05*(D-r) 

E<==  B 
F<==  10 

(Osing  the  value  of  PI ,  the  horizontal  increment 
plotting  and  labeling  is  determined.) 

If  PI  =  512  then 
3  <  =  64 

Else  If  PI  =  256  then 
G<  —  32 

Else  If  PI  =  64  then 
G<  =  8 

Else 

3<  —  4 
End  If 
Plot  B,C,1 

(This  lines  and  places  tic  merits  on  the  horizontal  ax 
Do  for  I  *  E  to  0  by  -G 
3  Plot  I*  (B-A| /E+A,C,2 

Plot  I*  (B-AJ  /E+A,0  (D-D /150r2 
S-;  Plot  I*  (B-  A)  /  E+  A, C  r  2 


End  Do 


(This  lines  and  places  tic  marts  on  tha  vertical  a xrs.) 
Do  for  I  =  0  to  ? 

Plot  A,I*(D-C)/F+C,2 
Plot  A ♦  (B-A)/150#r*(D-C) /?♦;,  2 
Plot  A  ,  I*  ( D-C)  /?*•  2,2 
End  Do 

C  har  act  er_  Size  1.2,1,. 7,0 
(This  numerically  libels  tha  varricai  axis.) 

Do  for  I  =  F  to  0  by  -1 

Plot  A-.075*  (B-A)  ,1*  (D-T)  /F+C,1 


Label  I/F 


End  Do 


(This  sets  up  a  character  string  for  labeling  -he  vertical 


axis. ) 


If  P3  =  0  then 


AS  TRI  NG<==  'DATA  POINTS'' 


L<  =  1 


Else  If  P3  *  2  than 


ASTRING<«  ' L INES/MIIRD RADIAN ' 


L<*»  1  /  (2*L) 


End  If 


(This  numerically  libels  tha  horizontal  axis.) 


Do  for  I  *  0  to  E  by  3 


V.  CONCLUSION 


The  work  reports!  in  this  thesis  supports  the  model 
predicted  and  measured  in  Crittenden,  and  others,  for  the 
long  exposure  case  [Bef.  2].  The  vidicon,  in  replacing  the 
mechanical  slit  scanning  system,  shows  no  degrading  of  the 
signal  data.  It  also  demonstrates  that  a  good  approximation 
for  the  point  spread  function  may  be  made  by  recording  a 
single  TV  line  through  a  laser  spot.  This  line  is  then  used 
to  calculate  the  one-dimensional  line  spread  function.  The 
linearity  of  the  video  tape  recorder  is  seen  in  the  results 
of  the  MTF  and  Cn2.  The  463  oscilloscope  is  the  workhorse 
for  the  entire  systea.  It  effectively  displays,  stores, 
digitizes,  and  transfers  data  on  a  real  tiae  basis. 

The  overall  system  does  not  have  the  capability  of  the 
system  described  by  Crager.  However,  the  relative 
simplicity  cf  the  structured  programming  technique  coupled 
with  a  lower  equipment  cost  demonstrates  that  comparable 
measurements  and  data  evaluation  using  this  equipment  can  be 
made.  Further  investigation  should  include  use  of  the  Data 
Precision  6000  digital  waveform  analyzer  thereby  allowing 
data  to  be  sampled  and  analyzed  on  a  near  real  time 
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APPENDIX  A 


0 
1 
2 

3 

4 

5 

6 

7 

8 
9 
10: 
11: 
12: 
13: 
14: 
15: 
16: 
17: 
18: 
19: 
20: 
21: 
22: 
23: 
24: 
25: 
26: 
27: 
28: 
29: 
.30: 
.21: 
.32: 
33: 
34: 
35: 
.26 : 
37: 
38: 


"HE-NE: SOURCE  ON  9, SCALE  ON  ?,&  COLLIMATED  ON  5"  : 
"GA-AS: SOURCE  ON  8 , SCALE  ON  6  .&COLLIMATED  ON  4": 
ent  "SELECT  CODE  FOR  PRINTER" ,A;dev  "print"  ,A 
fmt  l,z  ,c  ;tmt  :2,f  8  .1  ,z;  trot  3  ,£8.  1 
dim  18(10241  ,R(512]  , Y [ 5121  ,C8(20]  ( 20 ]  ,E8(20] 

aim  1(512]  ,  AS  [  32  J  ,B$  [2  0J  ,A  (  2,2]  ,B (2,1]  ,C[  2,1] 
but  "DATA"  ,I8,.3  ?ina  Y  ,1  ,R;0*R+W+V+U ;  1+L;buf  "I",  1 8, 3 
beep; asp  "ENTER  LASER  USED";wait  1500 
ent  "HE-NE=1 ,GA-AS=0"  ,W 
it  W=l;  "HE-NE "+E8, 9+0,  3^p  2 
"GA-AS "+E 8; 8+0 

beep;dsp  E8&"  SOURCE  DATA  RECORDED?" ; wai t  1500 
ent  "1=YES ,Q=NO"  fV 

if  V=0  and  n=l;  sfg  G;cll  "  'TRANSFER  '  (0)  ;Dmp  2 
if  V.=0  and  w=u;£fg  0  jell  '  'TRANSFER  '  ( Q ) 
beep;ent  "PLOT  OF  SOURCE?, 1=YES  ,0=NC"  ,R 
if  R*l;Cll  'PLOT '( 512 ,Q  ,  G) 

beep; asp  E8&"  SCALE  DATA  RECORDED?" ;wait  15C0 
oeep;ent  "1=YES , G*NO" , U * 

ir  U=0  ana  w=l;sfg  0;cll  'TRANSFER '  (7+0)  ;jmp  2 
if  U=G  ana  w=G; sfg  0 ;  ell  ’  'TRANSFER  '( 6+0) 
beep ;asp  "WHAT  IS  SCALE  FACTOR?" ; wai t  1000 
ent  "SCALE  FACTOR*? ;  0  GETS  PLOT"  ,S 

U#0;u-2+wcll  'PLOT  '(512, Or  0) 

U=0 ;  ell  'PLOT '(512, 0,0) 

-3 


ir  S=0  ana 
if  S=G  ana 
if  a=0;gmp 

S+L 

if  W=l;  9*.T*0;^nip 
8+.T+0 

laf  T,I$,R[*J 
pos ( 18 , "R: " ) +B 
val (1$ (B+2 ,B+4] ) 
pos  ( 1 8 ,  "  I T : ". )  +B 
I8(B+.3,£+4]*A8 
laf  T-2 , 1 8 ,R ( * ] 
pcs(I8,"R;")+B 
val(I8(B+2,B+4]  )*r2 
pos ( 18 , " IT : " ) +B 
I  $ (B+3 ,B+4 ] +B8 


rl 
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39 

40 

41 

42 

43 

44 

45 

46 

47 
4b 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 


laf  .T-2,I$,R{*] 
pos ( 1$ , "R: " ) ♦£ 
val(I$(B+2,B+4] )+r3 
pos(I$,  "LT:'\)*B 
IS[B+3,B+4]  +C$ 

if  r2  =  rl  and  A$.=B$  ana  r2=r3  ana  B$=C$;jmp  2 

beep;asp  "ERROR  IN  TIKE  SCALE “;stp 

HB":ld£  Q,  1$  ,R[*  ] ;  ina  Y;ara  R+Y ;  Q+A+B+C 

for  1*1  tc  512 

A+Y  [  I  ]  -*A ;  if  I  <51;B+Y[  I]  *B 

if  I  >463  ;B+Y[*I  ]  ♦B 

next  I?prt  A ,B 

B/100*D;A-512*D-A 

for  1*1  tc  512;Y{I]  -D+i { I ] ;  C+Y  ( I  ]  -*C 
if  C<A/2;  I-*J 
next  I;prt  C,J 
for  K=1  tc  256 

Yf  ( J.-K+l)ir.od512+l]  +  Y [  (K+J-l)  mcc512  +  l]  ♦I[K] 
next  K;ina  R;ara  I*R 

"IMAGE  POINT  SPREAD  FCN"-I$;rcr  1Q,I$,R[*] 

beep;ent  "PSF  CU.TPUT  =  1  AND/OR  CON!"  ,Z 

beep ; ent  "PSF  PLOT=l  AND/OR  CONT"  ,Y 

if  Z*l;cll  'TABLE '(10) 

if  Y=l;cli  "PLOT' (22, 1C, 1) 

ir.in  (I  [*]  )  ♦r22;max  (I  [  *].)  ♦r23 

prt  r22,r23 

oeep;ent  "LSF  OUTPUT.*  1  AND/CR  CCNT"  ,2 
beep; ent  "LSF  FLOT=l  AND/CR  CONT"  ,Y 

Cll  "LSF  "IMAGE  LINE  SPREAD  FCN”-I$;rcf  10,I$,R[*] 
if  Z*l;cll  ‘'.TABLE* '(10) 
if  Y  =  l;cll  "PLOT '(32 ,10,1) 

"CALCULATION  CF  FXFORM  OF  LSF";laf  10 , 1$ ,R(*J ; ina  I 
ara  R+I 

L(1J*R[1] ; for  1*2  tc  256 
I [ I ] *R[I]*R(514-I] ;next  I 
R  [  256]  *R[257]  ;rcf  10,I5fR[*] 

Cll  "FXFORM "FXFORM  OF  IMAGE  LSF " +1  $ ; rc f  10,I$,R[*] 
beep;ent  "FXFORM  OF  LSF  OUT’PUT.*l  AND/CR  CONT"  ,Z 
beep;ent  "FXFORM  OF  LSF  PLOT*l  AND/OR  CCNT”  ,Y 
if  Z*l; cll  "TABLE" (10) 
if  Y*l;cll  "PLOT"  (32,10,2) 

"B*D I A  OF  CBSCUR/DIA  OF  OBJ  LENS  IN  METERS"; 

"B*0 . 064/. 164 " ; 

"rl*SCALE  OF  DATA  IN  MICRORADIANS , rl  +  1. 60" ; 

"O-DIA  OF  OBJ  LENS  IN  METERS ?O*0 . 1 64M" ; 

"W -WAVELENGTH  IN  METERS"; 
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"CALC  CF  DIFFRACTION  LIMIT  POINT  SPREAD  FCN  : 
ina  R;it  W=l; 6.328e-7*W; jmp  2 
9 .  0  5e-7*K 

. 064/. 164+B; . 164*0  ;L*rl 

B  ~2*D  ;  1-D*H  ?  7T*rl*0*le- 6/Vv*Z 

tor  L=1  tc  256 

Z* ( 1-1) *Y ; it  Y>3G;gtc  +3 

("AIRY  '  (Y)  -D*  'AIRY"  (Y*B)  )  *2/H''2-R[I] 

I*K;gto  +2 

0  *R [ I ] 

txa  .5;asp  Rll];next  I;"OPTICS  FUNCTTCN"*I  $ 
rci  11 ,1$ ,R[*J 

beep; ent  “OPTICS  FCN  PLOT=l  AND/OR  CONT"  ,Y 
it  Y=l;cll  "PLOT"  (32 ,11,1) 

Stg  5;cll  "LSF";  "LSF  OF  OPTICS  FCN  "*I$;  rci  11,I$,R[*J 
beep  ;ent  "PLOT  OF  LSF  OF  OPTICS  FCN  =  1  AND/OR  CONT" ,Y 
xt  Y=i  ;c  11  "  'PLOT"  (  32 , 11,  1; 
sfg  5;ina  I;ara  R*I ;  I.{  1] +R[  1] 
for  1=2  to  256;I  [ I ] *R [ I ] *R [514-1] ;next  I 
R[25b] *R(257J ;rcr  il,I$,R(*J 

cl 1  "FXFCRM"; "FXFCRM  CF  LSF  OF  OPTICS" *16; tct  ll,I$,R[*j 
beep  ;ent  "PLOT  OF  FXFM  OF  GPTICS  =  1  AND/OR  CONT" ,  Y 
if  Y*l;Cll  "PLOT '(256,11,2) 

"CALCULATION  OF  QUOTIENT  OF  TWO  FOURIER  TRANSFORMS": 

laf  10 ,1  $ ,R[* ] ; ina  I;ara  R*I 

Idf  11 , 1  $ ,R ( * ] 

for  1=1  tc  256 

I [I ]/R(I] *R[I] ;next  I 

"  M'F  CF  SYSTEM"  *18;  rc  t  10,I$,R[*J 

beep  ;ent  "PLOT  OF  MTF  OF  SYS=1  AND/CR  CONT",Y 

it  Y-l;cll  "PLOT '(32,10,1) 

0*rl7*r 18*rl9*r  20*r21*r22*r2  3*r24*r2  5 

tor  1=1  to  96 

(  ( 1-1)  *L )  ~  ( 10/3)  *rl8 

rl8+rl7*r 17 

( ( I — 1 ) *L) * (5/3) *r 19 

rl9+r20*r  20 

I*r2 1 

if  R(L]  <=0;  gtc  "CC" 

(In (R(I] ) *r24) *rl9+r 23*r2 3 
r 24+r 25*r25 
"CC ":next  I 

prt  rl7,r20,r21,r23  ,r25 

r20*A( 1,1] ;r  21*A( 1 ,2] ; rl7*A[2  ,1]  ;r20*A[2,2] 
r25*C (1,1]  ;  r23*C. [2,1] 
inv  A*A;mat  A*C*B;flt  5;aprt  B 
"CALCULATION  OF  CNSC" : 

B[l,l]/(-21.49*145*r20*Sv"(.-.33333)  )  *r22 
prt  "CNSQ*"  ,r22;  stp 


134:  cfg  ?laf  2 

135:  "IRANSFER" : beep; asp  "ENSURE  EQUIP.  SET  UP  PROPERLY " 

136:  wait  1500 

137:  asp  "PRESS  CONTINUE  WHEN  READY" ; stp 
138:  dsp  "waiting" ;wait  1500 
139:  but  "DATA" ;  0-*Z 
140:  tfr  703, "DATA" 

141:  dsp  "transferring" 

142:  rds { "DATA" ) *E; if  E=-l;jrcp  0 
143:  clr  703 

144:  dsp  "setting  output" 

145:  pos(I$,"%")-A;l+N 

146:  for  I=A+2  tc  687  ty  16 

147:  for  J=1  tc  16;if  U  +  I  >688 ;  0-*>r  J  ;  jmp  3 

148:  num  (I$[I+J] )  ♦rJ 

149:  rJ-*-Y(N] 

150:  N+1*N 
151:  next  J 
152:  next  I 

153:  min(Y[*J  )  +E;for  1  =  1  tc  512;  Y  (I  ]  -£-*Y  [  I  ] ;  next  I 
154:  ina  R;ara  Y-»R;  rcr  pl,I$  ,R [ * ] ; ret 
155:  "LSF"  :  if  tlg5;laf  11 ,1$  ,R(  *j  ;  jrop  2 
156:  ldf  10 , 1  $ ,R [* ] 

157:  ina  I;ara  R-*I;ror  1=1  tc  24 
158:  1+J ;  asp  I ;  I  £  I  ]  -C. 

159:  y(I*I+J*J)-R 

160:  2*  (  (1-trc  (R)  )  *  I  ( int(F)  ]+frc(R)*I[int(R)+l])+g-*-Q 
161:  J  +  l->J;if  R<24;jir.p  -2 
162:  Q-L[I] 

163:  next  I ; ina  R;ara  I+R; if  tlg5;rcf  11, I$,R[*] ;crg  ;ret  pi 
164:  ret  pi 

165:  "FXFORM" :raa;9+N; ina  I ; if  rlg5;laf  11 , 1$ , R ( * ] ; jmp  2 
166:  laf  1Q,I5,R[*] 

167:  tt/2"  (N-l)  *T 

168:  for  M=1  tc  N 2*  (M-ft)  ♦rO 

169:  for  J=0  to  2  *  (M-l).-l;  ell  '  'BI  '  ( J  ,P  ,N-1) 

170:  cos  (P*T)-*C;sin(P*T)*(l-2*.tlg7)+P 
171:  for  I=2*rO*J+l  tc  2*r0*J+rG 
172:  R(I]-rl;R(I+rO]-r2 
173:  I[IJ-r3;I£I+rOJ*r4 

174:  rl+r  2*C+r  4*P-*-R(  I]  ;  r  3+r4*C-r2*P*I  [  I] 

175:  rl-r2*C-r4*P-*Rl  I  +  rC]  ;  r  3-r  4*C  +  r2*  P- 1  [  I  +  rQ  ] 

176:  next  I;next  J;asp  M;next  M 
177:  for  1=0  to  2"‘N-l?cll  'fll  '( I  ,J  ,N ) 

178:  if  I-J>0;gto  "BE" 

179:  if  I=J ;gto  "INC" 

180:  R(I+1]/V(2"N)-P;I(I  +  1]/V  (2"N)-Z 
181:  R[J+1]*R[I+1J  ;I  £  J+l]  -  I  £1+1] 

182:  P+R { J+l] ;Z+I(J+1] 
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183 

184 

185 
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204 
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207 
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211 
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214 
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"  I N  C  "  :  R  ( I  + 1  ]  /  v  ( 2  )  ♦R  { I  +1  ]  ;I  [1  +  1]  /  v  <  2  *1  [  1  +  1 ! 

"B£" : nex  t  I 

it  tig  5; ret  11 , 1  $ , R [ * ]  ; a eg ; re t 
deg ;  r  et 

"£I " : 0  *p2 ; p l+p4 
tor  Z=1  tc  p3 
p4/2-*-p4 ;2*p2*p2 
it  tre (p4) #0; p2+l-»-p2 
int(p4) *p4 
next  Z ;  ret 

"AIRY":  it  pl<C;beep;d£p  "errcr-argument<0" ;stp 

it  pl=0;  l-*-r4;  r  e  t  r4 

G-*-r5;it  pl>15;3mp  2 

2G+10*pl-pl-'2/3-*-r6 ;  jir.p  =2 

90+pl/2-»r6 

it  pl<  5;  6  +  pl-*-rl2;  jrr.p  2 

1.4*pl+6C/pl-*-rl2 

rr.ax(int  <rl2)  ,int(3+pi/4)  )+rl2 

tor  &=rl2  to  r6  by  3 ;le-2o ♦ro;  C*rl2*r  14 

stg  10;  it  iY/2=int(t:/  2)  ;ctg  10 

for  J=1  tc  M-2 ; 2* (K-J ) *r 8/pl- rl3+rl 5 ; rb*r 13 

rl5-*-ra,*ir  M-J-2=C;  rl5-*-r4 

cmf  10;  rl4+2*  r8*t lglG+r 14;  next  0 

2*r8/pl-rl3-rl5 

rl4+r!5-*r  14;  r4/rl4*r4 

it  abs  (r  4-rS).-aos  (r  4*le-6)  <=0;  2* r4/pl+r4 ; r et  r4 
r4+r5;next  h 

beep; dsp  "ACCURACY  NOT  OBTAINED  "  ;  v,e i  t  1500  ;ret  r< 

"TABLE"  :iat  p],I8,R[*J 

tor  1=1  tc  32;rcr  J=1  tc  15 

wrt  "print.  2"  ,R  1 16  ( I -1) +0  ] ;  nex  t  J 

wrt  "pr  int.  3"  ,R  [  lbl]  ;  next  I;ret  0-*Z 

"PLOT" : lat  p2,I$,R[*] 

0-*-A;pl*o;nu.n(  R  l  *  ]  )  *C ;  .T,a  x  ( R  [  *  ]  )  *D 

SCI  A-.l(B-A)  ,B+> .  C  5  «B-A)  ,C-.  l(L-C)  ,0-4.0  5  (D-C) 

B-*-E ;  10-*-F 
it  pl=512;64-G 
if  pl=256;32+G 
it  pl=64;8-*-G 
if  pl=32;4-*G 
pit  B,C,1 

tor  I=E  tc  0  by  -G 

pit  I  ( B-A)  /E+A  ,C  ,  2 

pit  I  (B-A) /E+A  ,C+ (0-0/150,2 

pit  I (B-A) /E+A ,C, 2 

next  I 

tor  1=0  tc  F 

pit  A, I ( D-C) /F+C  ,  2 

pit  A+(B-A)/150,I(D-C) /F+C  ,2 
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pit  A, I (C-C) /F+C,  2 
next  I; pen 
csiz  1.2,1,.  .7 
f  xa  1 

for  I=F  tc  0  by  -1 

pit  A-.075 (B-A) , I (D-C)/F+C ,1 

lbl  I/F 

next  I 

fxc  0 

it  p3=Q; "DATA  POINTS" ♦A?; 1*L 

if  ,p3  =  2;  "LINES/MICRGRADI  AN" -*-A$ ;  1/(2*L)  ♦L 

for  1=0  tc  E  by  G 

pit  A+(  I/E- .  025.)  (B-A)  ,C-.025(D-C)  ,1 
lbl  I*L 

next  I;  if  p3#0  and  p3f  2;  "MICRCRADIANS"-*-AS 
if  p2  =  2;S-L 

pit  .4  (B-A)  +A  0  5(D-C)  +C  ;  iol  AS 

pit  0  7  (B-A) +A  ,  ..3  (D-C ) +C;  csiz  1.2,1, .7,9 

"NORMALIZED  INTENSITY" ♦A?;  lbl  AS 

csiz  .  5 , 1 , 1. 5 , 0;  0  +  1 

pit  I ,max ( R ( * ] ) ,1 

for  1=1  tc  pi ; pl t  1-1 ,R[I ] ; next  I 

csiz  1.2,1,..  7,0 

if  p2#Q;gtc  +2 

beep ;ent  "PLOT  LABEL?", IS 

pit  . 6B  , .  9D  ,  1;  lbl  "PLOT  OF" 

pit  .6B, .670,1; lbl  IS 

pit  .  bB  ,  .  8  4C  ,  1 ;  i  tl  ES,"  LASER" 

pen;ctg  ;ret  U*Y 
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MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS-  1963-A 


0:  "COMPUTED  SOURCE  PATTERN" : 

Is  ina  R;  10C0*A  ;t2*L-C ;  trk  ljfaf  10 

2:  tor  1=1  tc  256 

3:  L*<L-])*F 

4:  F*2/ ( 2*C*2)  *G 

5:  it  G>13;jmp  2 

6:  A*exp(.-G)  ♦Rf  I]  ;I*K;  jmp  2 

7:  0 ♦R ( I ] 

8:  next  I 


10 

11 

12 

12 

14 

15 

16 

17 

18 

19 

20 
21 
22 


'"COMPUTED  SOURCE"-*T $ 
ref  10,I$,R(*J 

beepjent  "PLOT  OF  COMPUTED  SOURCE=l  AND/OR  CCNT",f 
it  y=l;cll  "'PLOT ' (32,10,0) 

"CONVERTS  PSF  10  LSF " : 

:  Cll  'LSF  ' 

!  “LSF  OF  COMPUTED  SOURCE"*I$ 
ref  10 , 1$ ,R [ *] 

:  beepjent  "PLOT  Or  LSF  OF  COMP  SC0RCE*1  AND/OR  CONI" 

:  if  y=ljcll  "PLOT '(32, 10,0) 

"CALCULATION  OF  FXFORM  OF  LSF  OF  COMPUTED  SOURCE": 

ina  I;ara  R-*I 

K1HI1] 

tor  1=2  tc  256 


.if 
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24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


1. 1 1  ] -*R  ( I  ]  ♦R  [  514-1  ]  ;next  I 
R(256]-Rf'257]  ;rct  10,16, Rl*] 

Cll  "FXFORM' 

“FXFORM  OF  LSF  OF  CALC  SOURCE” -1 5 ; ret  1G,I$,RI*J 
beep ;ent  "PLOT  OF  FXFM  OF  CALC  SRCE=1  AND/CR  CGNT" ,  Y 
if  Y=ljCll  "PLOT'(256,10,2) 

"CALCULATION  OF  PRODUCT  OF  2  FX FORMS  ”  : 

laf  10 , 1  $  ,  R  l  *]  ;  ina  I;ara  R+IjtrK  0 

lot  10,I$,R[*] 

for  1=1  to  256 

R [ 1 1  *L[I1-R(  I]  j  next  I 

ref  10 , 1$ ,R[ *] 

"TRANSFER  FCN  OF  SOURCE*ATMCSPHERE*OP.TICS "  : 
lot  11 , 1  $ ,  R  ( *  J j  ina  I;ara  R+I 
Ill]*R(l}-R[l] 
for  1*2  to  256 


35 

40 

41 

42 

43 

44 

45 


R [ I ] * I ( I ] ♦R 1 1 ] -R [ 514-1 ] 
next  I 

R [ 256] ♦Rf  257] 

"FXFM  OF  SGURCE*ATMOS*OPTICS"*I$ 
ref  10 , 1$ ,R( *] 

beepjent  "PLOT  OF  FXFM  PRODUCTS*l  AND/OR  CONT"  ,Y 
if  y»l;cll  "PLOT '(32,10,2) 
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"INVERSE  FXFCRM  GIVES  TARGET  LSF  : 

Stg  7;cil  ' 'FXFCRM erg  7;"INVERSE  FXFCRM I $ 

ref  10,I$,R[*j 

ceep ;ent  "PLOT  OF  INV  FXFM=i  AND/GR  CCNT"  ,1 
if  Y=l;cll  ‘'PLOT '(32,10,0) 

"CONVERTS  LSF  TO  PSF  BY  ABEL  TRANSFORM": 

Cll  'ABEL';  "ABEL  XFORM  PSF"-I$ 
ref  10 , 1$ ,Rl *] ; sf g  0 

beep;ent  "PLOT  OF  ABEL  XFGRM=1  AND/CR  CCNT" ,Y 
if  Y=1  ;cll  * 'PLOT '(32 ,10,  C)  ;stp 

"CALCULATES  FRACTION  OF  POWER  INSIDE  CIRCLE  OF  RADIUS  R 

.  25*  7r^R  £1  ]  -R  J 1 J 

for  1=2  to  256 

2n*I*R ( I ] +R( L-l] ♦R [I  ] 

next  I 

"POWER  FRCTN  INSIDE  CIRCLE" +1  $ 
ref  lU,I$,R[*j 

beep; ent  "PLOT  OF  POWER  =  l  AND/CR  CCNT"  ,Y 

if  Y=l;cll  "PLOT  '  (32  ,10,0) 

stp 

"LSF"  :  ir  rig  £;  lui  11 ,1$  ,R(*J  ;o*P  2 
laf  10 , 1$ ,R [ *] 

ina  I;ara  R*I;for  1=1  tc  24 
l*J;asp  I ; I ( I J ♦O 
y (I*I+J*J)*R 

2*  ( (l-rrc(R) )  *1  (int(R)  ]  +rrc  (R)  *i.  [  int  <R)  +1]  )-Kl*Q 
0  + 1*0 ;  if  R<24;iic.p  -2 

u*im 

next  l;ina  R;ara  I+R; if  rlg5;rcr  ) 1,1 $,R[*] ;crg  ;ret  pi 
ret  pi 

"FXFCRM"  :  raa;S-*N;  ina  I ;  if  fig  5;  lot  11 , 1  $  ,R  [  *  ]  ; ^mp  2 
lot  lu , IS  »R( *] 
ir/2'(N-l)  +T 

for  M=1  to  N; 2" (N-M) *r0 
for  J=0  tc  2*  (M-l).-l ;  cll  " '31  '  ( J  ,P  ,N-1) 
.ccs(P*T)*C;sin(P*T)  *  (l-2*flg7)>P 
for  I=2*rG*J+l  tc  2*r0*J+r0 
R[I]-rl;R(I+r0] -r2 
I [IJ >r3; L [I+r0] +r4 

rl+r2*C+r4*P+R(I] ; r 3+r4*C- r2*P*I 1 1 ] 

rl-r2*C-r4*P*R{  I+r0]  ;  r3-r4*C  +  r  2*P->I  [I+r0] 

next  I;  next  J;asp  M;next  M 

for  1=0  tc  2-N-l;cll  'fil'd, J,N) 

if  L-0  >0;  gto  "BB" 

if  I=J;gto  "INC" 

R(I  +  1]  /y  (2"N)*P;L(I  +  l]/y(2''N)+Z 
RtJ+1]  -RH  +  1]  ;I  (J+1J  -I [1  +  13 
P*R  [  J+1J  ;Z+I(J41J 

"INC"  :R[I+1)  /V  (2*N)  *R[I+1)  ;I  H  +  l)  /y  (2~N)  ♦![  I+1J 


95:  "EB":next  I 

96:  cegjret 

57:  "BI" : Q*p2j pl*p4 

98:  for  2*1  to  p3 

99:  p4/2-*p4;  2*p2-*p2 

100:  if  frc(p4) #0;p2+l*p2 

101:  int(p4)-*>p4 

102:  next  Z;ret 

103:  "ABEL"  :  ina  R;ldf  10  ,1$  ,R[*]  ;  R[l] +N 
104:  1.4*R(l)-1.8*R[2]*h.4*R[3]*R(l] 

105:  for  1*2  to  64 
106:  R{IJ*M 

107:  ,4*N+.2*M-.6*R[I+1] *R{I] ;M*N;next  I 

108:  for  1*1  to  64 

109:  Rfl  ]/  (2*V  C  (!■*»•  1)  “2-I*L) ) -R(I  ] 

110:  for  J*I+1  to  64 

111:  RUJ+RUl/V  (  (JH..1)  '2-IH)*R[I] 

112:  next  J 

313:  if  R[L]<.01;0-R[I] 

114:  R  [  I  ]  /ir*R  (  L]  ;dsp  Ijnext  I 

115:  for  1=65  to  512; G*R [I J ; next  I; ret 

116:  "PLOT” : lat  p2,I$,R[*] 

117:  0-A;pl*B;min(R[*] ) +C ; max (R [*] ) 

118:  scl  A  - .  1  (B-A)  ,B-h.G5  «B-A)  ,C-.1(D-C)  ,D+.05(D-C) 

119:  B*E; 10+F; if  flgOjprt  A,B,C,D 

120:  it  pl=512;  64*G 

121:  if  pl=256;32*G 

122:  if  pl*64;8*G 

123:  if  pl=32;4*G 

124:  pit  B,C,1 

125:  for  l=E  to  0  by  -G 

126:  pit  I (B-A) /E+A  ,C , 2 

127:  pit  I  (B-A)  /E+A  ,C+  (D.-C)/ 150 , 2 

128:  pit  I (B-A) /E+A, C, 2 

129:  next  I 

130:  for  1*0  to  F 

131:  pit  A , I ( D-C ) / F  +C , 2 

132:  pit  A+(B-A)/15Q ,1 (D-C)/F+C ,2 

133:  pit  A, I ( L-C) /F+C , 2 

134:  next  I;pen 

135:  csix  1.2,1, .7 

136:  fxa  1 

137:  for  I*F  to  0  by  -1 

138:  pit  A-.075(&-A) ,1 (D-C)/F+C,l 

139:  lbl  I/F 

140:  next  1 

141:  fxa  0 

142:  if  p3*2; "LINES/MICRCRADIAN"*A$; l/( 2*L)  +L 
143:  for  1*0  to  E  by  G 
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144:  pit  A+(I/E-.C2S} (B-A) ,C-.025(D-C) ,1 
145:  lbl  I*L 

146:  next  I; it  p3#2; "MICRORADIANS" *A$ 

147:  it  p3«2;S-*L 

148:  pit  .4(B-A)+A^-.05(D-C)+C?lbl  A$ 

149:  pit  .-.07(B-A)+A,..3(D-C)+C;csiz  1.2,1,.7,9C 
150:  "NORMALIZED  INTENSITY" +A$; lbl  A$ 

151:  csiz  .  5, 1,1.5, 0;Q-*T 

152:  pit  I ,max(R[*] ) ,1 

153:  for  I»1  to  pl;plt  I-1,R[I J ; nex t  I 

154:  csiz  1.2,1,. 7,0 

155:  pit  .  6£  , .  9D ,  1;  lol  "PLOT  OF" 

156:  pit  .6B,.87D,1; lbl  1$ 

157:  pit  .6B , .84D, 1; lbl  E$ LASER" 

158:  "RET" :pen;ctg  ;ret  0+Y 
*17524 
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